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Liquid  water  transport  in  gas  channels  may  influence  and  even  prevent  gas  supply  in  a  proton  exchange 
membrane  (PEM)  fuel  cell,  so  significantly  affect  cell  operation.  A  two-phase  two-dimensional  lattice 
Boltzmann  model  is  employed  in  this  paper  to  simulate  the  development  and  interaction  of  two  liq¬ 
uid  droplets  in  a  gas  channel  of  a  PEM  fuel  cell,  focusing  mainly  on  parametric  effects,  including  the 
gas  flow  velocity,  initial  droplet  distance,  different  micropore  combinations,  and  the  gas  diffusion  layer 
(GDL)  surface  wetting  properties,  on  liquid  droplet  transport  behaviors.  Results  confirm  that  an  increased 
gas  stream  velocity  and  liquid  pore  distance  may  prevent  liquid  droplet  interaction  and  enhance  liq¬ 
uid  water  removal.  Numerical  simulations  further  indicate  that  different  pore  size  combinations  may 
promote  droplet  interaction,  particularly  with  a  large  pore  in  front  of  a  small  one.  On  the  contrary,  a 
more  hydrophobic  GDL  surface  can  decrease  liquid  droplet  interaction,  benefit  liquid  water  removal,  and 
consequently  improve  PEM  fuel  cell  water  management. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Liquid  water  transport  and  distribution  is  a  paramount  issue 
in  proton  exchange  membrane  (PEM)  fuel  cell  water  management 
and  plays  a  dictating  role  in  fuel  cell  performance  and  durability.  In 
the  past  decade,  many  experimental  [1-10]  and  numerical  studies 
[11-31]  have  been  conducted  in  this  area  to  elucidate  the  funda¬ 
mental  liquid  water  transport  mechanisms  and  their  effects  on  fuel 
cell  operation.  In  recent  years,  to  obtain  deeper  insight  into  the 
underlying  physics,  the  lattice  Boltzmann  method  has  also  been 
employed  for  direct  numerical  simulations  of  liquid  water  trans¬ 
port  in  porous  materials  of  PEM  fuel  cells  with  realistic  micropore 
structures  [32-35].  These  meso-scale  studies  paint  a  clearer  picture 
of  liquid  water  transport  process  in  PEM  fuel  cells  with  physical 
details  not  attainable  using  current  experimental  techniques. 

Although  good  progress  has  been  achieved  in  numerical  mod¬ 
eling  of  liquid  water  transport  in  the  porous  materials  in  PEM 
fuel  cells,  there  still  remains  a  weak  area  concerning  liquid  water 
dynamic  behaviors  in  gas  channels.  Liquid  water  transport  in  gas 
channels  may  influence  and  even  prevent  gas  supply  in  a  PEM  fuel 
cell  and  thus  significantly  affect  fuel  cell  operation.  No  comprehen¬ 
sive  and  robust  numerical  model  capable  of  accurately  handling 
liquid  water  transport  in  the  gas  channels  of  a  PEM  fuel  cell  exists 
in  the  open  literature.  So  far,  almost  all  the  available  PEM  fuel  cell 
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models  treat  liquid  water  transport  in  a  gas  channel  as  a  mist  flow, 
which  assumes  tiny  liquid  droplets  moving  synchronously  with  the 
gas  stream.  The  mist  flow  model  can  only  cover  a  very  small  portion 
of  the  liquid  water  transport  regions  in  the  gas  channels  of  PEM  fuel 
cells. 

Given  the  practical  importance  of  liquid  water  transport  behav¬ 
iors  in  PEM  fuel  cell  gas  channels,  many  experimental  studies  have 
been  carried  out  to  lay  a  foundation  for  fundamental  understanding 
of  the  phenomena  [1-3].  Yang  et  al.  built  an  optically  transparent 
PEM  fuel  cell  to  observe  liquid  water  transport  processes  in  the 
cathode  gas  channel  [2].  It  was  found  that  liquid  water  droplets 
could  emerge  from  the  gas  diffusion  layer  (GDL)  interior  to  its  chan¬ 
nel  surface  at  preferential  locations,  with  droplet  diameters  ranging 
from  a  few  to  hundreds  of  micrometers.  These  liquid  droplets  would 
subsequently  grow  and  interact  on  the  GDL  surface,  and  under  cer¬ 
tain  operation  conditions,  could  block  gas  flow  and  cause  possible 
cell  operation  shut-down.  Zhang  et  al.  [3]  further  studied  liquid 
water  removal  phenomena  in  the  cathode  gas  channel  using  a  sim¬ 
ilar  transparent  fuel  cell.  The  experimental  results  revealed  two 
major  modes  of  liquid  water  removal  from  the  GDL  surface;  one 
being  droplet  detachment  by  the  shear  forces  of  gas  flow,  which 
subsequently  led  to  a  mist  flow  in  the  gas  channel,  and  the  other 
being  capillary  wicking  of  liquid  droplets  onto  the  hydrophilic 
channel  walls,  which  resulted  in  the  annular  film  flow  and/or  liquid 
slug  flow. 

To  obtain  more  physical  insights,  a  number  of  numerical 
simulations  have  been  further  performed  to  study  liquid  water 
transport  mechanisms  in  PEM  fuel  cell  gas  channels.  Quan  et  al. 
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Nomenclature 

c  lattice  space  and  time  ratio 

cs  lattice  speed  of  sound 

e  lattice  velocity  vector 

/  density  distribution  function 

F  external  force 

Fads  adhesive  force 

^int  inter-particle  force 

Fg  body  force 

G  constant  controlling  interaction  strength 

P  pressure 

A P  pressure  difference,  Pa 

r  droplet  radius,  m 

R  gas  constant 

s  index  function 

T  temperature 

t  time 

u  macroscopic  fluid  velocity  vector 

ug  inlet  gas  stream  velocity,  m  s-1 
x  space  position  vector 

Greek  symbols 
6  contact  angle 

p  fluid  density 

p0  constant  density 

pw  tunable  parameter 

a  surface  tension,  N  m-1 

r  relaxation  time 

ooa  weighting  factor 

if/  potential  function 

\//o  constant 

Superscripts 
eq  equilibrium 

Subscripts 

a  direction  index 

g  gas  phase 

int  inter-particle  parameter 


[36]  investigated  liquid  water  droplet  behaviors  and  interaction 
with  the  air  stream  and  channel  walls  in  a  serpentine  gas  channel 
of  a  PEM  fuel  cell.  The  volume  of  fluid  (VOF)  method  was  applied 
to  track  the  liquid  water  surface.  Numerical  results  indicated  that 
the  bend  area  of  the  serpentine  channel  would  exert  significant 
influences  on  the  air  flow  and  liquid  water  transport  and  distribu¬ 
tion  inside  the  channel.  It  was  found  that  under  specific  operation 
conditions,  liquid  water  flooding  might  block  the  gas  channel 
and  prevent  reactant  supply.  The  same  method  has  later  been 
employed  by  this  research  group  to  study  liquid  water  transport  in 
straight  parallel  channels  with  the  inlet  and  outlet  manifolds  in  a 
PEM  fuel  cell  stack  [37].  Recently,  Hao  and  Cheng  [38]  used  the  lat¬ 
tice  Boltzmann  method  to  simulate  liquid  water  droplet  dynamic 
behaviors  on  a  hydrophobic  GDL  surface.  The  detailed  process  of 
a  liquid  water  droplet  emerging  from  a  micropore  in  the  GDL,  its 
continuous  growth  on  the  GDL  surface,  and  the  final  detachment 
from  the  micropore  was  clearly  presented.  It  was  concluded  that 
the  gas  flow  velocity  and  GDL  surface  wetting  properties  would 
strongly  affect  the  liquid  droplet  movement.  This  research  work 
also  demonstrated  that  the  lattice  Boltzmann  method  could  be  a 
powerful  tool  for  numerical  simulations  of  liquid  water  transport 
phenomena  in  the  gas  channels  of  PEM  fuel  cells. 


Fig.  1.  Schematic  of  the  D2Q9  scheme,  showing  fixed  lattices  in  a  2D  space  and  9 
specified  particle  streaming  directions  (including  a  stationary  one). 


In  this  paper,  a  two-phase  two-dimensional  lattice  Boltzmann 
model  is  employed  for  direct  numerical  simulations  of  the  emer¬ 
gence,  growth,  and  subsequent  interaction  of  two  liquid  water 
droplets  on  the  GDL  surface  in  a  gas  channel  of  a  PEM  fuel  cell. 
Particular  attention  has  been  paid  to  the  effects  of  a  number  of  key 
influential  parameters  on  liquid  water  droplets  development  and 
interaction,  including  the  gas  flow  velocity,  initial  droplet  distance, 
different  micropore  combinations,  and  the  GDL  surface  wetting 
properties.  The  simulation  results  may  help  to  gain  more  funda¬ 
mental  insights  on  liquid  water  dynamic  behaviors  in  PEM  fuel  cell 
gas  channels. 

2.  Theoretical  formulation 

In  this  paper,  we  implemented  a  lattice  Boltzmann  model  (LBM) 
in  a  two-dimensional  configuration  to  simulate  liquid  droplets 
development  and  interaction  in  a  gas  channel  of  a  PEM  fuel  cell.  LBM 
is  a  powerful  and  convenient  tool  for  simulating  multi-phase  flows, 
since  it  can  maintain  phase  interface  automatically  during  numer¬ 
ical  calculation  without  resorting  to  more  complicated  interface 
tracking  techniques  commonly  used  in  the  conventional  numerical 
methods  for  multi-phase  simulations,  as  in  the  VOF  approach.  The 
lattice  Boltzmann  method  simulates  fluid  flows  through  the  evolu¬ 
tion  of  the  particle  density  distribution  functions  on  fixed  lattices. 
The  calculations  proceed  in  two  steps:  collision  and  streaming. 
The  collision  procedure  relaxes  the  density  distribution  functions 
on  a  lattice  toward  the  equilibrium  states,  while  the  streaming 
step  moves  the  particle  distributions  to  the  immediate  neighboring 
lattices  with  specified  velocity  magnitudes  and  directions.  A  num¬ 
ber  of  different  multi-phase  LBM  approaches  have  been  proposed 
and  have  been  successfully  applied  in  many  fluid  flow  simulations 
[39-44],  including  liquid  water  transport  in  the  porous  materials 
of  PEM  fuel  cells  [32-35]. 

The  present  numerical  calculations  are  based  on  the  lattice 
Boltzmann  method  with  a  single  relaxation  time  collision  oper¬ 
ator  (the  so-called  LBGK  model)  and  the  popular  D2Q9  scheme, 
as  shown  in  Fig.  1.  The  Shan-Chen  multi-phase  model  is  further 
implemented  for  studying  two-phase  liquid  droplets  development 
and  interaction  [39,40,43].  The  lattice  Boltzmann  equation  can  be 
expressed  as 

fa(x  +  eaAt,  t  +  At)=/„(X,  t)-  (1) 

where /is  the  density  distribution  function,  x  refers  to  space  posi¬ 
tion,  t  time,  a  the  velocity  direction,  and  r  the  relaxation  time.  It 
should  be  emphasized  that  the  lattice  Boltzmann  model  is  formu¬ 
lated  and  solved  in  dimensionless  form,  so  the  relevant  parameters 
in  this  model  are  all  dimensionless  ones  without  units  [39,40,44]. 
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In  Eq.  (1),  the  equilibrium  distribution  function  is  calculated  as 


C(X,  t)  =  coap 


1  + 


ea  •  uei 


(e~g  -u^)2 
2  cf 


(uiq)2 

2Cj 


(2) 


More  details  concerning  this  two-phase  lattice  Boltzmann 
model  can  be  found  in  the  open  literature  [39-44]. 

3.  Model  validation 


In  Eq.  (2),p  =  ^2  afa  is  the  fluid  density.  In  the  D2Q9  scheme,  the 
weighting  factor,  wa,  is  chosen  to  be  4/9  for  a  =  0, 1/9  for  a  =  (1,  2, 
3, 4),  and  1  /36  for  a  =  (5,  6,  7,  8)  to  fully  recover  the  Navier-Stokes 
equation  [44].  The  parameter,  ea,  represents  the  9  discrete  veloci¬ 
ties,  as  shown  in  Fig.  1,  and  are  defined  as 

[e0,  e\ ,  e2,  e3,  e4,  e5,  e6,  e7,  e8] 

_  [oio— 1  0  1—1—1  l] 

—  c  [o  0  1  0  —11  1  —1  — 1 J  w 

where  c=  Ax/ At  is  the  ratio  between  the  lattice  space  distance  and 
time  step.  The  parameter,  cs  in  Eq.  (2),  is  given  as  c/V 3  in  the  LBGK 
model  [44].  In  this  dimensionless  LBGK  simulation,  Ax  and  At  are 
both  set  to  be  unity. 

In  the  Shan-Chen  model  [39,40],  extra  forces  have  been 
introduced  to  account  for  fluid  inter-particle  and  fluid-solid  inter¬ 
actions.  These  forces  are  incorporated  into  the  equilibrium  velocity 


as 

tF 

ueci  =  u'  +  — 
p 

(4) 

where  the  macroscopic  velocity  is 

u,  _  ^cfaea 

P 

(5) 

The  extra  forces  in  Eq.  (4)  is  defined  as 

F  —  ^int  +  Fads  +  Fg 

(6) 

In  Eq.  (6),  the  parameter,  Fg,  represents  a  uniform  steady  body 
force.  The  inter-particle  cohesive  force  is  introduced  as  follows: 


Fint(X>  0  =  ~G\KX,  +  6a  At,  t)ea  (7) 

a 

where  G  is  a  constant  representing  the  interaction  strength 
between  neighboring  particles,  while  the  weighting  factor  coa  has 
the  same  values  as  in  Eq.  (2).  The  parameter,  i/q  is  a  density- 
dependent  potential  function  given  as 

Vf(/°)  =  Vfoexp(A^°)  (8) 

where  the  parameters,  t/T)  and  p0,  are  two  constants. 

In  this  study,  the  following  equation  of  state,  which  allows  for 
liquid-gas  phase  coexistence,  is  applied: 

P=pRT+^x/r2(p)  (9) 

where  the  dimensionless  parameter,  RT ,  is  equal  to  1/3,  according 
to  the  gas  kinetic  theory  [41  ]. 

The  parameters,  G,  t/T),  and  p0,  are  chosen  to  be  -125,  4,  and 
200,  respectively,  based  on  previous  two-phase  studies  using  the 
LBGK  model  [40,41  ]  and  our  requirement  for  large  gas  and  liquid 
density  ratio  in  the  present  simulation. 

In  Eq.  (6),  an  adsorption  force  is  also  included  to  account  for 
liquid  and  solid  phase  interactions,  which  is  formulated  as 

Fads (X,  t)  =  -G\j/{p{x,  t))^<w«Vf(Av)s(x  +  ea  At,  t)ea  (10) 

a 

where  s  is  an  index  function  with  a  value  of  1  if  the  site  at  x  +  eaAt 
is  in  solid  phase,  with  a  value  of  0  otherwise.  pw  is  a  free  parameter 
for  tuning  different  wall  properties,  e.g.  the  contact  angle. 


The  two-phase  lattice  Boltzmann  model  briefly  described  in  the 
last  section  has  been  programmed  into  a  computational  software 
package  in  our  group  and  validated  herein  using  two  test  cases.  The 
first  test  case  concerns  liquid-gas  interaction  and  phase  separa¬ 
tion,  which  would  lead  to  surface  tension.  Fig.  2  shows  many  liquid 
droplets  with  different  diameters  obtained  after  conducting  phase 
separation  calculation.  The  results  provide  necessary  information 
for  determining  the  surface  tension.  According  to  the  Laplace  equa¬ 
tion,  the  following  formula  exists: 

A  P=-  (11) 

r 

which  basically  states  that  the  pressure  difference,  A P,  inside  and 
outside  a  droplet  is  proportional  to  the  surface  tension,  cr,  and 
inversely  proportional  to  the  droplet  radius,  r. 

Using  the  results  in  Fig.  2,  the  droplet  surface  tension  can  be 
calculated  using  pressure  differences  and  droplet  radii,  and  the  cal¬ 
culated  results  are  presented  in  Fig.  3.  The  surface  tension  obtained 
from  the  present  LBM  approach  agrees  very  well  with  the  Laplace 
equation,  as  clearly  verified  in  Fig.  3. 

The  second  test  case  concerns  liquid-solid  interaction,  which 
would  lead  to  different  wetting  characteristics  of  a  solid  surface, 
e.g.  different  contact  angles  as  illustrated  in  Fig.  4.  By  changing  the 
parameter,  pw,  in  Eq.  (10),  the  strength  of  liquid  and  solid  interac¬ 
tions  can  be  varied  in  the  present  two-phase  LBM  approach.  This 
naturally  results  in  different  solid  surface  contact  angles,  as  illus¬ 
trated  in  Fig.  5,  which  shows  a  liquid  droplet  in  contact  with  a  solid 
wall  under  different  wetting  conditions. 

The  numerical  results  obtained  in  the  two  test  cases  clearly 
indicate  that  the  two-phase  lattice  Boltzmann  model  presented  in 
the  last  section  is  capable  of  handling  liquid-gas  and  liquid-solid 
interactions.  Therefore,  it  is  suitable  for  simulating  liquid  droplets 
development  and  interaction  in  a  gas  channel  in  a  PEM  fuel  cell. 

It  should  be  noted  that  because  of  the  lack  of  high-resolution 
experimental  data  in  both  time  and  space  concerning  the  fol¬ 
lowing  two  droplets  interaction  in  the  gas  channel  of  a  PEM  fuel 
cell,  model  validation  in  this  paper  is  conducted  against  two 
representative  two-phase  phenomena.  Although  this  is  a  common 
practice  and  has  been  proven  to  be  a  reliable  approach  [34,38-44], 
direct  validation  using  high-resolution  experimental  data  can 
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Fig.  2.  Liquid  droplets  with  different  diameters  after  a  phase  separation  calculation. 
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provide  further  confidence  and  should  be  carried  out  once  the 
experimental  data  become  available. 

4.  Results  and  discussion 

A  two-phase  two-dimensional  lattice  Boltzmann  model  has 
been  presented  and  validated  in  the  earlier  sections,  and  is 
employed  here  to  study  liquid  droplet  development  and  interac¬ 
tion  in  a  gas  channel  of  a  PEM  fuel  cell,  as  illustrated  in  Fig.  6. 
Previous  experimental  results  [2,3]  have  clearly  illuminated  that 
liquid  droplets  would  emerge,  move,  and  interact  on  the  GDL  sur¬ 
face  in  the  PEM  fuel  cell  gas  channel.  The  present  numerical  studies 
are  intended  to  provide  more  details  on  liquid  droplet  interaction 


processes  and  help  enhancing  fundamental  understanding  of  the 
dynamic  phenomena.  The  present  LBM  simulations  focus  mainly 
on  the  effects  of  various  key  parameters,  including  the  gas  stream 
velocity,  micropore  size  and  distance,  and  GDL  surface  wetting 
properties,  on  liquid  droplet  transport  behaviors. 

As  shown  in  Fig.  6,  the  gas  channel  depth  is  chosen  to  be  500  p.m, 
and  its  length  is  1000  p.m.  There  are  two  micropores  in  the  GDL 
on  one  side  of  the  channel,  with  the  pore  width  chosen  as  30  p.m 
unless  otherwise  specified.  These  are  representative  parameters  for 
a  typical  PEM  fuel  cell.  In  the  simulations,  liquid  water  is  forced  into 
the  gas  channel  through  the  two  micropores  with  a  specified  inlet 
velocity  of  0.08 ms-1  [38].  Gas  stream  flows  into  the  gas  channel 
with  a  parabolic  shape,  assuming  a  fully  developed  laminar  flow. 
The  GDL  surface  contact  angle  is  specified  at  120°.  Based  on  grid 
independence  studies,  500  x  250  regular  lattices  are  placed  in  the 
gas  channel  region  in  the  x  and  y  directions,  respectively,  while 
15  x  50  lattices  are  used  in  each  of  the  micropores. 

In  the  following  LBM  simulations,  effects  of  the  inlet  gas  stream 
velocity  on  droplet  interactions  are  first  investigated.  In  this  part  of 
simulations,  the  distance  between  the  two  micropores  are  specified 
at  240  fxm,  while  the  averaged  inlet  gas  stream  velocity  varies  from 
2.0  to  4.0ms_1. 

The  transient  variations  of  the  two  liquid  droplets  transport  are 
displayed  in  Fig.  7.  The  movement  of  the  two  droplets  is  mainly 
influenced  by  the  shear  force  from  the  gas  stream,  droplet  surface 
tension,  and  two  droplets  interaction.  As  shown  in  Fig.  7,  with  a  low 
inlet  gas  velocity  at  2.0  m  s-1 ,  the  two  droplets  emerging  from  the 
micropores  would  be  slightly  elongated  along  the  flow  direction, 
but  neither  of  them  could  be  pulled  away  from  the  pores  under  the 
relatively  weak  shear  drag.  Eventually,  the  leading  droplet  would 
grow  sufficiently  large  to  merge  with  the  trailing  droplet,  forming 
a  single  large  deforming  droplet  under  the  influence  of  the  shear 
force.  The  existence  of  this  large  droplet  would  further  affect  the 
gas  stream  velocity  distribution,  as  shown  in  Fig.  8. 

As  the  inlet  gas  velocity  is  increased  to  3.0  m  s-1 ,  the  two  emerg¬ 
ing  droplets  would  exhibit  different  interaction  behaviors.  As  the 
shear  drag  acting  on  the  leading  droplet  becomes  stronger,  this 
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Fig.  5.  Liquid  droplet  attached  to  a  solid  wall  at  different  contact  angles. 
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droplet  is  detached  from  its  feeding  pore  and  subsequently  collides 
with  the  trailing  droplet.  The  liquid  water  region  formed  by  the  two 
droplets  exhibits  a  quite  flat  shape  under  the  strong  gas  stream 
shear  force.  As  shown  in  Fig.  7,  once  the  inlet  gas  stream  velocity  is 
further  increased  to  4.0ms_1,  the  droplet  merging  phenomenon 


would  disappear.  Under  a  large  shear  drag,  the  leading  droplet 
would  be  pulled  away  from  its  micropore  prior  to  3.44  ms,  but  it 
could  not  catch  the  trailing  droplet,  which  is  also  detached  from  its 
feeding  source  very  quickly.  Afterwards,  the  two  droplets  would 
move  separately  on  the  GDL  surface  in  the  gas  channel.  Results  in 


t  =  2.20  ms  t  =  2.40  ms  t  =  2.50  ms 


t  =  4.12ms  t  =  1.60ms  t  =  1.20ms 


t  =  5.18ms  t  =  4.04ms  t  =  1.44ms 


t  =  5.8  ms 


t  =  4.14ms 


t  =  1.64ms 


ug  =  2.0  ms  1 


ug  =  1.0ms  1 


ug  =  4.0  ms  1 


Fig.  7.  Effects  of  inlet  gas  stream  velocity  on  liquid  droplet  interaction. 
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Fig.  7  confirm  that  a  high  gas  stream  velocity  is  helpful  to  liquid 
droplet  removal  in  a  gas  channel  of  a  PEM  fuel  cell,  which  would  in 
turn  benefit  the  gas  stream  flows,  as  illustrated  in  Fig.  8. 

Effects  of  the  two  micropores  distance  on  liquid  droplet  inter¬ 
action  are  next  studied.  In  these  three  cases,  the  inlet  gas  stream 
velocity  is  specified  at  4.0  m  s-1 ,  while  the  pore  distances  are  set  at 


200, 240,  and  280  pan,  respectively.  Detailed  simulation  results  are 
presented  in  Fig.  9. 

At  a  relatively  short  distance  of  200  pan,  the  leading  droplet 
would  be  detached  from  its  feeding  pore  under  a  strong  shear  drag 
from  the  gas  stream,  but  comparing  with  the  other  two  cases,  its 
detaching  time  is  slightly  delayed  by  the  close  influence  of  the 
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Fig.  9.  Effects  of  micro  pore  distance  on  liquid  droplet  interaction. 
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Fig.  10.  Effects  of  different  micro  pore  combinations  on  liquid  droplet  interaction. 


second  droplet.  The  leading  droplet  would  exert  strong  influence 
on  the  trailing  droplet  by  reducing  the  shear  force  it  confronts.  In 
this  way,  the  leading  droplet  could  catch  up  to  and  merge  with  the 
trailing  droplet  before  latter  is  pulled  away  from  the  micropore. 

As  the  pore  distance  is  increased  to  240  p.m,  the  effect  of  the 
leading  droplet  on  the  trailing  droplet  becomes  greatly  weakened. 
As  discussed  earlier  in  Fig.  7,  the  trailing  droplet,  in  this  case,  could 
be  pulled  away  from  its  micropore  soon  after  the  leading  droplet  is 
detached.  The  two  droplets  would  then  move  separately  in  the  gas 
channel. 

When  the  pore  distance  is  further  increased  to  280  pun,  inter¬ 
action  between  the  two  droplets  essentially  disappears.  As  shown 
in  Fig.  9,  the  two  droplets  would  be  detached  at  almost  the  same 
instant  and  behave  independently  in  the  entire  process.  This  sit¬ 
uation  should  be  beneficial  for  liquid  water  removal  and  water 
management  in  a  PEM  fuel  cell. 

The  pore  size  effects  on  liquid  droplets  interaction  in  a  gas 
channel  are  also  investigated  in  this  paper.  In  the  simulations,  the 
inlet  gas  velocity  is  4.0ms_1,  the  pore  distance  is  set  at  240  pan, 
but  the  two  micropores  are  given  different  sizes.  As  shown  in 
Fig.  10,  one  of  the  micropore  is  reduced  to  20  pan,  leading  to 


different  micropore  combinations  as  20-20,  20-30,  30-20,  and  30- 
30.  Results  from  the  last  30-30  case  have  already  been  presented 
and  discussed  in  Figs.  7  and  9. 

Results  from  the  20-20  case  are  similar  to  those  in  the  30-30 
case.  The  two  droplets  are  initially  elongated  slightly  along  the  flow 
direction.  The  leading  droplet  would  then  be  pulled  away  from  its 
feeding  pore,  but  before  it  could  collide  with  the  trailing  droplet, 
the  second  droplet  would  also  become  detached.  During  the  entire 
simulation,  the  two  droplets  move  separately  in  the  gas  channel. 

For  the  20-30  case  in  Fig.  10,  the  existence  of  the  leading  small 
droplet  alters  the  shear  force  direction  acting  on  the  trailing  large 
droplet,  and  consequently  the  trailing  droplet  is  strongly  stretched 
along  the  flow  direction.  Confronting  a  strong  incoming  flow,  the 
small  leading  droplet  could  be  easily  detached  from  its  microp¬ 
ore  before  the  trailing  large  droplet  could  be  pulled  away.  The  two 
droplets  would  be  able  to  merge  thereafter. 

With  a  large  pore  in  front  of  a  small  one,  droplet  interaction 
would  behave  much  differently.  In  the  30-20  case,  the  leading  large 
droplet  is  more  stretched  than  the  trailing  small  droplet.  The  large 
and  much-elongated  leading  droplet  could  easily  touch  the  small 
trailing  droplet  and  merge  with  it  very  quickly.  During  this  process, 
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Fig.  11.  Effects  of  surface  contact  angle  on  liquid  droplet  interaction. 


however,  both  droplets  cannot  be  pulled  away  from  their  feeding 
pores,  and  the  merged  liquid  droplet  becomes  strongly  stretched 
along  the  flow  direction. 

Results  in  Fig.  10  clearly  show  that  different  pore  size  com¬ 
binations  could  exert  significant  influences  on  liquid  droplet 
development  in  a  gas  channel.  Placing  two  micropores  with  differ¬ 
ent  diameters  together  would  generally  lead  to  enhanced  droplet 
interaction,  particularly  with  a  large  pore  in  front  of  a  small  one. 

Effects  of  the  GDL  surface  wetting  property  on  liquid  droplets 
development  and  interaction  are  finally  simulated.  In  these  cases, 
the  inlet  gas  stream  velocity  is  set  at  2.0  m  s-1 ,  while  the  pore  dis¬ 
tance  is  240  fxm. 


The  effect  of  a  hydrophilic  GDL  surface  with  a  contact  angle  of 
30°  on  droplet  interaction  is  illustrated  in  Fig.  1 1 .  With  a  hydrophilic 
surface,  the  two  liquid  water  regions  emerging  from  the  micropores 
are  closely  attached  to  the  channel  wall  as  thin  films,  spread¬ 
ing  along  the  flow  direction.  As  a  result,  the  two  liquid  regions 
could  merge  together  very  quickly,  comparing  with  the  hydropho¬ 
bic  cases  in  Figs.  7  and  11.  The  merged  large  liquid  region  remains 
attached  to  the  feeding  pores  as  a  flat  liquid  film  and  becomes  very 
difficult  to  eject. 

For  a  hydrophobic  GDL  surface,  liquid  droplets  behave  dif¬ 
ferently  with  an  increased  GDL  surface  contact  angle  at  150°, 
comparing  with  the  results  in  Fig.  7  with  a  contact  angle  at  120°. 
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As  the  surface  contact  angle  increases,  the  interaction  process 
between  the  two  droplets  becomes  weaker,  which  is  attributable  to 
the  interaction  between  the  liquid  droplet  and  solid  wall.  As  such, 
liquid  droplet  removal  may  be  enhanced  by  a  more  hydrophobic 
GDL  surface. 

Results  in  Figs.  7  and  1 1  indicate  that  a  strongly  hydrophobic 
GDL  surface  may  decrease  liquid  droplet  interaction,  benefit  liquid 
water  removal,  and  consequently,  improve  water  management  in 
a  PEM  fuel  cell. 

5.  Conclusions 

In  this  paper,  the  Shan-Chen  two-phase  lattice  Boltzmann 
model,  along  with  the  popular  D2Q9  scheme,  has  been  employed  to 
simulate  two  liquid  droplets  development  and  interaction  in  a  gas 
channel  of  a  PEM  fuel  cell.  The  two-phase  two-dimensional  lattice 
Boltzmann  model  has  been  validated  using  two  representative  test 
cases  concerning  phase  separation  and  liquid-solid  interaction.  In 
present  numerical  simulations,  particular  attention  has  been  paid 
to  the  effects  of  a  number  of  key  influential  parameters  on  liquid 
water  droplet  transport  behaviors,  including  the  gas  flow  velocity, 
initial  droplet  distance,  different  micropore  combinations,  and  the 
GDL  surface  wetting  properties. 

Results  confirm  that  an  increased  gas  stream  velocity  and  liq¬ 
uid  pore  distance  may  prevent  liquid  droplet  interaction,  enhance 
liquid  droplet  removal,  and  benefit  water  management  in  the  gas 
channel  of  a  PEM  fuel  cell.  Numerical  simulations  further  indi¬ 
cate  that  different  pore  size  combinations  may  promote  droplet 
interaction,  particularly  with  a  large  pore  in  front  of  a  small  one. 
On  the  contrary,  a  more  hydrophobic  GDL  surface  can  decrease 
liquid  droplet  interaction,  benefit  liquid  water  removal,  and  con¬ 
sequently  improve  PEM  fuel  cell  water  management.  The  present 
lattice  Boltzmann  simulations  are  intended  to  provide  more  details 
on  liquid  droplet  interaction  behaviors  and  enhance  fundamental 
understanding  of  the  related  dynamic  phenomena. 
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